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ABSTRACT 

A transiting extrasolar planet sequentially blocks off the light coming from the 
different parts of the disk of the host star in a time dependent manner. Because of 
the spin of the star, this produces an asymmetric distortion in the line profiles of 
the stellar spectrum, leading to an apparent anomaly in the radial velocity curves, 
known as the Rossiter - McLaughlin effect. Here, we derive approximate but 
accurate analytic formulae for the anomaly in the radial velocity curves, taking 
into account the stellar limb darkening. The formulae are particularly useful in 
extracting information on the projected angle between the planetary orbit axis 
and the stellar spin axis, A, and the projected stellar spin velocity, V sml s . We 
create mock samples for the radial curves for the transiting extrasolar system HD 
209458 and demonstrate that constraints on the spin parameters (V sin I s , A) can 
be significantly improved by combining our analytic template formulae and the 
precision velocity curves from high-resolution spectroscopic observations with 
8-10 m class telescopes. Thus, future observational exploration of transiting 
systems using the Rossiter - McLaughlin effect will be one of the most important 
probes for a better understanding of the origin of extrasolar planetary systems, 
especially the origin of their angular momentum. 

Subject headings: planets and satellites: individual (HD 209458b), techniques: 
spectroscopic 
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1. Introduction 

With more than 130 extrasolar planets discovered so far, major scientific purposes in this 
field are rapidly moving from mere detection to characterization of the planetary systems, i.e., 
statistics of planetary masses, orbital periods, eccentricities, stellar metallicities, and so on. 
Most of these systems have been discovered through the periodic change of radial velocities 
of the central stars. One of them, HD 209458b, was first discovered spectroscopically and 
was soon later found to exhibit a transit signature in front of the stellar disk for a duration 
of ~ 2 hr in its orbital period of 3.5 days (Henry et al. 2000; Charbonneau et al. 2000). More 
recently, a few additional extrasolar planets have been discovered photometrically from a 
transiting signature in their light curves and later confirmed spectroscopically (e.g., Udalski 
et al. 2002c; Bouchy et al. 2004; Alonso et al. 2004). 

Indeed, such transiting planets provide important information for the extrasolar plan- 
etary systems that is otherwise unavailable: planetary size, atmospheric composition, and 
the degree of (mis) alignment of the planetary orbit axis and the stellar spin axis. Among 
others, Queloz et al. (2000) showed that the planetary orbit and the stellar rotation of the 
HD 209458 system share the same direction; the planet sequentially blocks off the light from 
the approaching and then from the receding parts of the stellar surface. This produces a 
distortion in the line profiles of the stellar spectrum during the transit in a time-dependent 
manner, leading to an anomaly of the radial velocity curves, previously known as the Rossiter 
- McLaughlin (RM) effect in eclipsing binary stars (Rossiter 1924; McLaughlin 1924; Kopal 
1990). The signature of the RM effect was first reported by Schlesinger (1909) and the effect 
was later isolated and measured by Rossiter (1924) and McLaughlin (1924). The first 
attempt at a theoretical study of the RM effect was made by Petrie (1938) and the analysis 
was subsequently extended by Kopal (1942, 1945) to incorporate the effects of limb and 
gravity darkening and rotational and tidal distortion. In the context of extrasolar planetary 
systems, Queloz et al. (2000) numerically computed the expected amplitude of the radial 
velocity anomaly caused by the RM effect. They first put constraints on the stellar spin 
angular velocity and its direction angle with respect to the planetary orbit by comparing the 
observed velocity anomalies. Incidentally, Snellen (2004) proposed an interesting suggestion 
to use the RM effect as an efficient diagnostic of the atmosphere of a transiting planet. 

Obviously, this methodology provides unique and fundamental clues to understanding 
the formation process of extrasolar planetary systems. Planets are supposed to form in the 
proto planetary disk surrounding the proto star (e.g., Pollack et al. 1996). Thus, the stellar 
spin and the planetary orbital axes are expected to be aligned. In turn, any constraints 
on their (mis) alignment degree are useful clues to the origin of the angular momentum of 
planets and its subsequent evolution during possible migration of the planets into the close-in 
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orbits (e.g., Lin et al. 1996). In order to improve the reliability and precision of such results, 
analytic templates are of great value. In this paper, we show analytic formulae for radial 
velocities of transiting extrasolar planets. Those formulae can be used as standard templates 
in constraining a set of parameters for numerous transiting planets that will be detected in 
the near future with Corot and Kepler. 

The rest of the paper is organized as follows. Section 2 introduces a variety of parameters 
that characterize the dynamics of a planetary system and summarizes the radial velocity 
curve, neglecting the transit effect (or, equivalently, for a non rotating star). Section 3 shows 
a general theoretical framework of the RM effect for a star during the planetary transit. 
Sections 4 and 5 derive analytic expressions for radial velocity curves; in §4, we consider 
an idealistic case of the stellar intensity model without limb darkening and derive the exact 
analytic expressions. Notation used in the derivation of are summarized in Table 1. On this 
basis, §5 presents approximate formulae taking into account the stellar limb-darkening effect. 
We apply these analytic templates to a transiting system, HD 209458, in §6 and examine the 
sensitivity to the parameters of the system. Finally, §7 is devoted to the main conclusions 
and discussion. 



2. Radial velocity profile for a star with a nontransiting planet 

A close-in extrasolar planetary system may have multiple outer planets, but we focus 
here on a system that is well approximated by a two-body problem, i.e., one that consists of 
a central star (mass m s ) and a planet (mass m p ). Figure 1 shows the schematic configuration 
of the top view of the planetary orbit. The radial velocity curve of the star in the Kepler 
orbit can be described as follows (e.g., Murray & Dermott 1999). 

First, note that in the strictly two-body problem, the orbit of the planet with respect 
to the star is simply written as 

r = ° (1 - e2) (1) 
" 1 + ecos/' W 

where a is the semimajor axis, e is the eccentricity, and / is the true anomaly (angular 
coordinate measured from the pericenter direction). The true anomaly / is written in terms 
of the eccentric anomaly E, defined through the circumscribed circle that is concentric with 
the orbital ellipse as 

cosE — e 

cos/ = . 2 

J 1-ecosE v ; 
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If one introduces the mean motion n from the orbital period P orD of the system as 

n = -5—, (3) 

-'orb 

then E is related to the mean anomaly M as (Kepler's equation) 

M = E -esin£, (4) 

where M = n(t — r), where r is the time of pericenter passage. 

Using the parameters defined above, the radial velocity of the star along the line of sight 
of the observer (see Fig. 1) is written as (e.g., Murray & Dermott 1999) 

m„ nasini . . 

v ia d, s = ; 7= [sm(/ + tu) + e sin roj , (5) 

m s + m p VI - e 2 

where i denotes the inclination angle between the direction normal to the orbital plane and 
the observer's line of sight and we define — w as the longitude of the line of sight with respect 
to the pericenter (Fig.l). While / is not directly written as a function of the observer's 
time t, it is useful to rewrite equation (5) explicitly in terms of M — n(t — r) even in an 
approximate manner. For this purpose, one can use the following expansions with respect 
to the eccentricity e (e.g., Murray & Dermott 1999): 

00 1 d 

sin/ = 2 VI - e 2 V - — J k (ke) sin kM 

k de 

k=l 

= sinM + esin2M + e 2 Qsin3M- ^sinM^j +0(e 3 ), (6) 

2(1 _ e 2 ) ~ 
cos/ = — e H 2_\ ^ kike) cos kM 

e k=i 

9e 2 

= cosM + e(cos2M- 1) + — (cos 3M — cos M) + 0(e 3 ). (7) 

8 

Then equation (5) up to 0(e) should read 

m v na sin i r . . , , , . ^ r 

Vrad s ~ sm(M + w) + e sm(2M + w)\ . (8) 

m s + m p 



3. Radial velocity profile for a star with a transiting planet 



An occultation of a part of the rotating stellar surface during transit of the planet causes 
a time-dependent asymmetric feature in stellar emission/absorption line profiles. If the line 
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profile is not well resolved, the asymmetry results in an apparent shift of the central line 
position, which contributes additionally to the overall "observed" stellar radial velocity. In 
order to describe the effect quantitatively, we set the coordinate system centered at the star 
so that its y-axis is directed toward the observer (Fig. 2). The z-axis is chosen so that the 
stellar rotation axis lies on the y-z plane. We also define the angle A between the z-axis and 
the normal vector h p of the planetary orbit plane projected on the x-z plane (Fig.3), i.e., 

sin A sin i 

h p = [ cosi ] . (9) 
cos A sin % 

Then the position of the planet is given by 



v 



r p cos / 
.sin/ 


cos A sin(/ + zu) — sin A cos % cos(/ + zu) 
r p | sin i cos( f + zu) | , (10) 

sin A sin(/ + zu) — cos A cos % cos(/ + zu) 

where Rk(9) denotes the rotation matrix of an angle 6 around the fc-axis. 
In our configuration, the angular velocity of the star is given as 

fl s = (0, Q s cos I s , fl s sin I s ). (11) 

Then the velocity of a point R = (x, y, z) on the stellar surface is 

(z cos I s — y sin I s \ 
xsin/ s . (12) 

—X COS Is / 

Thus, the radiation at frequency v from that point suffers from the Doppler shift due to the 
stellar rotation by an amount 

Au fLrrsinie . , 
= — 13 

v c 

with respect to the observer located along the y-axis in the present case. 

Consider a specific (emission or absorption) line whose intensity at a point (x, z) on the 
projected stellar surface is given by I u (x,z) = I(x, z)H(u), where H(v) represents the line 
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profile. The observed flux is computed by integrating the Doppler-shifted intensity at each 
point over the entire (projected) surface of the star: 



F v = J + I{x,z)H{u-Au 



dxdz 
D 2 ' 



(14) 



where D is the distance between the star and the observer. The factor (1 + Au/u) 3 appears 
because of the Lorentz invariance of the quantity l v jv 3 . While our analysis is applicable to 
both emission and absorption lines, we consider an emission line centered at v — i/q in the 
following, just for definiteness. Then the line profile function satisfies 



/ 



H(v)dv = 1, (15) 
vH{y)dv = i/ . (16) 
Since H{y) is supposed to be sharply peaked only around v , we have approximately 



/ 



f{v)H(v)dvnf(v ), 



(17) 



for an arbitrary smooth function f{y). 



If the resolution of the observational spectrograph were sufficiently high, the line profiles 
of the star and the planetary shadow would be separated, or at least the asymmetric feature 
might be detected for transiting systems (e.g., Charbonneau et al. 1998, 1999). In reality, 
however, such a high spectral resolution is quite demanding, and here we assume a somewhat 
lower resolution. Thus, we simply compute the resulting time-dependent shift of the line- 
profile-weighted mean position v due to an asymmetric occultation of the stellar surface 
during the passage of the transiting planet. Using expression (14) and the properties of the 
line profile function (eqs. [15] to [17]), we obtain 



v = 



vF v dv 



J F v dv 



= "0 



1 + 



1 AuV Au T . NJ J 

1H l{x, z)dxdz 

vq J vq 

AzA 3 

1 H I I(x, z)dxdz 

u o J 



> . 



(18) 



Since the amplitude of the Doppler shift (eq. [13]) is small, one can safely expand v up to the 
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leading order of Ais/is as 



v = u 




(19) 



Therefore, the "apparent" radial stellar velocity anomaly due to the RM effect is expressed 
as 



Equation (20) is the basic relation between I(x, z) and Av s in our subsequent analysis. 
Figure 4 shows a schematic illustration of the RM effect. Depending on the inclination and 
the orbital rotation direction relative to the stellar spin axis, the velocity curve anomaly due 
to the RM effect exhibits rather different behavior. 

The remaining task is to evaluate the integrals adopting a certain model of a stellar 
surface intensity I(x,z). Note that previous literature in the analytic study of radial ve- 
locity curves focused on expressing the integrals in the radial velocity shift (20) in terms of 
Kopal's associated a-functions under some model assumptions of the stellar intensity (e.g., 
Hosokawa 1953; Kopal 1990). While such detailed and exact approaches are required for 
stellar eclipsing binaries, the evaluation of the a-function is a demanding numerical task. 
Furthermore, for those systems a variety of effects become important, including limb dark- 
ening, distortion of stars due to their rotation and tidal interaction, the reflection effect 
(heating by the radiant energy of the companion), and gravity darkening (variation of the 
surface brightness due to the local surface gravity acceleration change). For the extraso- 
lar planetary systems, on the other hand, the radius and mass of a planet are significantly 
smaller than those of the host star. Thus, most of those effects can be safely neglected, and 
one can derive simpler, and still practically useful, analytic formulae applying perturbative 
expansion. In what follows, we present such analytic expressions for the RM effect with and 
without the stellar limb darkening. 




(20) 



- 8- 



4. Analytic expressions for a uniform stellar disk (without limb darkening) 

As a step toward an analytic model for the RM effect for extrasolar planetary systems, 
let us consider first an idealistic case in which the limb-darkening effect is neglected. We also 
assume that the planet is completely optically thick and not rotating, which is also assumed 
in the next section. In this case, one can obtain the exact analytic expression even without 
the perturbative expansion. The intensity at (x, z) on the uniform stellar surface becomes 

Ho ;x 2 + z 2 < R 2 and (x - X p f + (z - Z p ) 2 > R\ 
I(x,z) = l , (21) 

t ; otherwise 

where X p = (X p , Y p , Z p ) is the position of the center of the planet and R s and R p denote the 
radii of the star and the planet, respectively. We evaluate equation (20) at complete transit, 
ingress, and egress phases in the following subsections (see Fig. 5). 



4.1. Complete transit phase 

During a complete transit phase, the position of the planet satisfies the relation [X 2 + Z 2 ^ ^ 2 < 
R s — R p . Thus, the range of the integral in equation (20) is simply given by the stellar surface 
area with the entire planetary disk sutracted, i.e., 

I I dxdz — > I dx I dz — I dx I dz. (22) 

J J J-Rs J-^/Rl-x 2 Jx p -R p J Z v -./Rl-(x-X v )' 2 



Then we obtain 



I(x,z)dxdz = n(R 2 s - R 2 p )I , (23) 
xl(x,z)dxdz = -X p ixR 2 p I Q . (24) 
Substituting these results into equation (20), we find 



Av s = n s X p sin Ly^-i ( 7 = f^) " (25) 

Equation (25) implies that the time dependence of the RM effect during the complete transit 
is entirely incorporated in the planet position, i.e., X p = X p (t). 
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4.2. Ingress and egress phases 

At ingress and egress phases, on the other hand, the location of the planet satisfies the 

1/2 

relation, R s — R p < [X 2 + < R s + R p . Just for computational convenience, we rotate 

the coordinates in a time-dependent manner so that the planet is always located along the 
new x-axis: 



X \ 1 / Xp Zp \ ( X 



z J R ay /X% + Z% V ~Zp X p I \ : 
Then the position of the planet is given by 



where 



(26) 



z p J V o 



(27) 



V P = ^±3 _ i (28) 

K s 

In the new coordinates, equation (21) is rewritten as 

io ; x 2 + z 2 < 1 and (5 - 1 - i] p ) 2 + z 2 > j 2 
I (.hi) { ' , (29) 

; otherwise 



and the moments of the intensity reduce to 



(30) 



J J I(x,z)dxdz = R 2 ^nl — JJ I(x, z)dzdx | , 

J J xl(x,z)dxdz = — ;j— p — J J (X p x — Z p z)I(x, z)dzdx, (31) 

where the range of the integrals denoted by S indicates the overlapping region between the 
stellar and the planetary disks and can be explicitly written as (dark shaded regions, Fig. 6): 



dzdx — > dx dz + dx ^^dz. (32) 

IS Jxo J-VT^W JX p -~f J-y/-y 2 -(x-X p ) 2 

Note that the planetary and stellar circles intersect at (x , ±z ), where 



■f-r% r , — -, \/('l 2 -vl)l(lp + 2) 2 -'< 2 ] 



x ° = l -WT^y " = ^-^ = ^r^) • (33) 
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Let us also introduce 

C = l + ^-^o= 2{1 + r]p) • (34) 

Physically speaking, this corresponds to the separation between the intersection and the 
center of the planet along the X axis, but we allow ( to be negative (see Fig. 6) as well. 
Then equations (30) and (31) are analytically integrated as 



I(x, z)dzdx = I Q [ sin 1 z - (1 + r] p )z + 7 2 cos l (C/l) ] , (35) 

S 



and 



J j (ZX P - zZ p )I(x, z)dzdx = I X p (l + rip) [-z ( + I 2 cos _1 (C/7)] , (36) 

respectively. 

Combining these results, we finally find that equation (20) reduces to 

Av s = Q s X p sin I s — — — — - 37) 

vr-sm 1 z + (1 + r]p) z o - 7 cos (C/7J 

While the above expression for Av s seems a bit complicated, this is the exact result for the 
radial velocity anomaly through relations (28), (33) and (34) in terms of the planet position 
(X p , Z p ) specified by equation (10). 



5. Effect of stellar limb darkening 

To be more realistic, we now take account of the effect of limb darkening, which pro- 
duces the radial dependence of the intensity of the stellar disk. Among the various models 
proposed so far (e.g., Claret 2000), we adopt a linear limb-darkening law as the simplest, 
but a practically realistic one. Introducing the linear limb-darkening coefficient e, the stellar 
intensity is now given by 

( J {1 - e(l - aO } ; x 2 + z 2 < R 2 and (x - X p ) 2 + (z - Z p ) 2 > R 2 p 
I(x,z) = l ' , (38) 

[ ; otherwise 

where \x is the cosine of the angle between the line of sight and the vector normal to the 
local stellar surface : 

" - f^W- (39) 
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With the limb-darkening effect, however, equation (20) can no longer be analytically inte- 
grated in an exact manner. Therefore, we construct approximate analytic formulae on the 
basis of the result without limb darkening (e = 0; see section 4). 



5.1. Complete transit phase 

Applying the analytic results of §4.1 to the stellar intensity model (38), equation (20) 
is formally rewritten as 

Av s = n s X p sinI s 7 2 1 n l 2 „ umv ( 4 °) 

where W\ and W 2 are defined as 

W 1 -- 



: JJ s dxdzy/l-(x* + Z*)/ff„ (41) 

W 2 = -^-^ JJ^dxdzx^l-(x 2 + z 2 )/Rl (42) 

The above integrals are carried out over the entire planetary disk. As discussed in Appendix 
A.l, they reduce to one-dimensional integrals, which can be expanded with respect to 7 = 
R p /R s . Specifically, equations (A10) and (A14) show perturbative expressions up to the 
fourth order in 7. The accuracy of the fourth-order perturbation expansion is within a few 
percent even for 7 ~ 0.3 (Fig. 15). In practice, however, the value of 7 is expected to be much 
smaller, 7 £ 0.1. In this case, higher order terms in equations (A10) and (A14) contribute 
merely ~ 1% to equation (40), and one can safely use 

W x {p) ~ 0, (43) 
W 2 (p) ~ (l-p 2 ) 1/2 , (44) 

where p = (X p 2 + Z 2 ) 1/2 / R s (0 < p < 1 - 7). 



5.2. Ingress and egress phases 

If the linear limb-darkening effect is taken into account, equation (37), describing the 
ingress and egress phases, now becomes 

Av s = Q S X P sin I s 
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(1 - e) {-z ( + 7 2 cos-^C/t)} + t - ~ ^4 
x l + Vp / 45 n 

7T (1 - |e) -(l-e){sin- 1 ^ -(l + ^)^o + 7 2 cos- 1 (C/7)}-e Ws' V ; 



where W 3 and W 4 are defined by 



W 3 = / / t&dSVl - x 2 - F, (46) 



W 4 = J J dxdz iVl-x 2 -z 2 . (47) 

Appendix A. 2 derives approximate analytic expressions (A23) and (A24) for equations 
(46) and (47), respectively, assuming that 7 1. Again, if 7 ;$ 0.1, they can be safely set 

as 

W 3 ~ 0, (48) 

^4 ^ I 7(7-0 ^ (Xc; ^ 7) W 2 (l- 7 ), (49) 
2 ^(1-7; -7, 7) 



where 



x c = x + ^y^, (50) 



/ \ f, • -1 / 7 2 -(^-l-^ P ) 2 l 

g(x;vp,v = (l-^)sm I ^— j> 



2 ^ 1/2 



+ ^{72 - (x - 1 - ^) 2 }{1 - x 2 - 7 2 + (z - 1 - Vp ) 2 }. (51) 

As shown in Figure 16, the accuracy of equations (48) and (49) is typically within a fractional 
error of 5%-10% percent. Nevertheless, their contribution to the total error budget for the 
velocity anomaly (45) is within a few percent (see §6.2). Thus, equations (48) and (49) are 
practically good approximations in most cases. 



6. Application to the HD 209458 system 

So far, HD 209458 is the only extrasolar planetary system in which the RM effect is 
detected; Queloz et al. (2000) reported the first detection of this effect with ELODIE spec- 
trograph on the 193 cm telescope of the Observatoire de Haute Provence. They numerically 
computed the radial velocity anomaly due to the RM effect for a variety of model parameters 
and compared these with the observed radial curves. They concluded that a = ±3?9 ^21° 
and V sin I s = 3.75 ± 1.25 km s" 1 , where a is the angle between the planet's orbital plane 
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and the star's apparent equatorial plane and V denotes the stellar surface velocity. These 
are written as a = cos _1 (sini cos A) and V = R S Q S according to the notation of our current 
paper. We summarize the current estimates of the stellar and planetary parameters for HD 
209458 in Table 2 and the best solution for the spin parameters by Queloz et al. (2000) in 
Table 3. Since HD 209458 remains the best target for the precise measurement of the RM 
effect, we consider in this section the extent to which one can improve the constraints on the 
spin parameters with our analytic formulae. 

6.1. Parameter dependence 

Adopting the linear limb-darkening law for the stellar intensity model, the RM effect for 
a system in the Keplerian orbit is specified by 10 parameters: the limb-darkening coefficient 
e, the orbital parameters of the system (a, e, i, w, and P or b), the size of the stellar and 
planetary disks (R s and R p ), the projected stellar surface velocity V sin I s , and the projected 
angle between the stellar spin axis and the normal direction of the orbital plane A. Except 
for the last two parameters (V sin I s and A) , these can be independently determined from the 
usual radial velocity and transiting photometric data, at least in principle. This is indeed 
the case for the HD 209458 system (see Table 2). Therefore, it is natural to ask about the 
extent to which one can put constraints on the two parameters V sin I s and A from the radial 
velocity anomaly during the transit due to the RM effect. 

Consider first the sensitivity to the spin parameters (Vsinl s , A). Figures 7 and 8 il- 
lustrate our approximations for Av s adopting the estimated parameters of the HD 209458 
system (Table 2) with and without the stellar limb darkening (i.e., e = and 0.64, specif- 
ically), respectively. The central transit epoch is chosen as t = 0. Then ingress starts at 
t = —1.55 hr, the complete transit lasts for —1.07 < t < 1.08 hr, and egress ends at t = 1.56 
hr for e = 0.1 (sometimes these four epochs are referred to as the first, second, third, and 
fourth contacts, respectively). The range of the spin parameters, Vsinl s and A, adopted in 
these figures roughly covers the uncertainties of the values of Queloz et al. (2000). 

Comparison of the two figures indicates that the radial velocity anomaly Av s is also 
sensitive to the linear limb-darkening coefficient e. Obviously the amplitude of the radial 
velocity shift Av r is sensitive to Vsinl s . The projected angle A shifts the zero point of 
the radial velocity anomaly at earlier (A > 0) and later (A < 0) epochs for the orbital 
inclination i < 90°. This produces an asymmetry of the shape of the radial velocity anomaly. 
Note that the behavior becomes opposite for the inclination i > 90°, corresponding to the 
parameter degeneracy between (i, A) and (180° — i, —A). Because of the different dependence 
of the overall radial velocity anomaly on the spin parameters (Vsinl s , A), one can put 
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more stringent constraints on those if our formulae are combined with future precision data 
attainable by 8 - 10 m class telescopes with a high dispersion spectrograph (HDS) such as 
Subaru HDS. 

Before addressing this issue in detail, it is helpful to clarify the dependence of the RM 
effect on the other remaining parameters. To investigate this, we quantify the variation of 
the radial velocity shift with respect to a specific parameter change p ~ p + dp by 



for p = a, R s , R p , P or b, V sin I s , e, and e. In practice, we systematically decrease the value of 
dp/p for each parameter down to 1CT 6 and ensure the convergence of the derivative. For the 
angular parameters p — i,zu, and A, we simply take their scaling values at 10°: 



Here, we confirm the convergence of the derivative by decreasing the value of dp/10° down to 
10~ 6 in these cases. Our analytic formulae are indeed useful in evaluating these quantities at 
the fiducial parameters of the HD 209458 system (Tables 2 and 3). The results are plotted 
in Figure 9 as a function of time. Note that our definition of SAv s is normalized by the 
fractional error in the parameter, i.e., dp/p. 

Figure 9 clearly shows that the stellar radius R s and the orbital parameters a, i, and 
P or b sensitively change the normalized radial velocity variation SAv s at the ingress and egress 
phases, while the spin parameters (V sin I s , X) have a relatively smaller effect on 5Av s . To 
quantify the actual deviation of the radial velocity shift caused by the systematic errors in 
the observation, we must multiply the observational uncertainty listed in Table 2. Then 
it turns out that the most sensitive parameter is Vsinl s , causing the 12m s -1 variation 
of the radial velocity shift. The other parameters change the radial velocity shift at the 
ingress and egress phases by less than 5m s _1 . This amplitude itself is comparable to the 
SAv s induced by the uncertainty of the projected angle A. Nevertheless, the different time- 
dependent effects of the various parameters can be used to break the parameter degeneracy, 
which would enable an accurate determination of the spin parameters A and Vsinl s . In 
addition, Figure 9 even suggests that a more precise determination of the orbital parameters 
other than the spin parameters is also possible by combining the RM effect with the usual 
radial velocity measurement. 



5Av s = lim 



Av s {p + dp) - Av s {p) 
dp/p 



(52) 



SAv s = lim 

o!p->0 



Av s (p + dp) - Av s (p) 
dp/10° 



(53) 
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6.2. Accuracy of our formulae 

While our analytic formulae presented in the previous section will improve the efficiency 
of the parameter estimations relative to fully numerical approaches, we have to address a 
couple of issues before applying them to real data: their accuracy and the effect of the 
finite exposure time. Our formulae with limb darkening are derived on the basis of an 
empirical approximation to the integrals of the stellar surface intensity (§5). Furthermore, 
the real data do not instantaneously sample the radial velocity, but are averaged over a finite 
exposure time. We directly test those effects against the numerical solutions of equation (20). 
Figures 10 and 11 compare three results: numerical integration of equation (20), our analytic 
formulae, and the numerical average of equation (20) over a At = 10 minute exposure time (in 
practice, we separately average the denominator and the numerator of the analytic formulae 
assuming At = 10 minutes and then take their ratio). These results are labeled A, B, and 
C and plotted in solid, dotted, and dashed curves, respectively, in the top panels. 

For e = 0, the analytic formulae (curve B) are exact, and the completely negligible 
difference A— B should be regarded as a welcome check of the accuracy of our numerical in- 
tegration scheme. The bottom panels suggest that the three results agree within an accuracy 
of ~ lm s _1 , which is below the typical radial velocity sensitivity achieved (~ 3m s _1 ) and 
is only comparable to the latest achievement by HARPS (Santos et al. 2004). Therefore, as 
far as the HD 209458 system is concerned, we can safely use our analytic formulae as useful 
templates for the RM effect even if a finite exposure time of the order of 10 minutes is taken 
into account. 



6.3. Mock analysis of the spin parameter estimation 

Now we are in a position to ask whether our formulae combined with precision spec- 
troscopic data can improve the previous constraints on the spin parameters (f2 s sin/ s , A), 
among others. For this purpose, we create mock data for the radial velocity anomaly of the 
HD 209458 system and fit them to the analytic formulae. Basically, the mock data were 
created adopting the central values of the parameters listed in Tables 2 and 3, but assigned 
an overall Gaussian random error of the rms amplitude 5 m s -1 , which is the level of preci- 
sion achieved with the Subaru HDS assuming an exposure time of At = 10 minutes (Winn 
et al. 2004). In light of the most recent sensitivity achieved by HARPS (Santos et al. 2004; 
~ lm s _1 ), our error assignment may still be conservative if the error is not dominated by 
other possible systematic errors. To mimic the effect of the finite exposure time, we numeri- 
cally integrate the denominator and the numerator of equation (20) separately over At = 10 
minutes. Then we take the ratio and assign the random error, as mentioned above. 
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Note that the number of independent data points during the transit phase (~ 3 hr 
including the ingress and egress phases) is 16 for At = 10 minutes. The generated mock 
data are then fitted to the analytic radial velocity anomaly to estimate the spin parameters. 
Here, the fitting is performed assuming prior knowledge of the remaining eight parameters. 

First, let us see how the spin parameters are reliably estimated from the x 2 fit- To 
examine this, we create 50, 000 mock realizations and calculate the joint probability distri- 
bution of the best-fit parameters under a certain prior knowledge of R s , R p , and a. We use 
the x 2 function, 



with a Av = 5m s _1 and N = 16. In this analysis, according to the result in Figure 9, we 
particularly focus on the five parameters Vsinl s , A, R s , Rp, and a. Their input values are 
3750 m s _1 , 0°, 1.146R©, 1.347Rj, and 0.0468 AU, respectively. We assume a set of different 
prior values for R s , R p , and a indicated in each panel of Figure 12 and then perform the x 2 
minimization over V sin I s and A. 

The results are plotted as contour levels in Figure 12. Here, the two contour curves in 
each panel represent the 1 and 2 a levels of the joint probability. The solid curves plotted 
along the horizontal and vertical axes are the probability distributions of Vsinl s and A, 
respectively, each of which is the projection of the joint probability distribution over the 
other parameter. Note that in both cases the probability distribution is well approximated 
by a Gaussian distribution with 1 a values of cr VsinJs ~ 310 m s _1 and a x ^ 3?4. This 
result indicates that the estimated value for V sin I s is rather sensitive to the assumed value 
of the planetary radius R p , while A can be estimated reliably even if R p is not known 
accurately. This comes from the fact that the velocity shift is roughly proportional to 
V sin I s Rp] however, the projected angle A is sensitive only to the asymmetry of the radial 
velocity shift curve. Another important aspect is that the uncertainty of the prior knowledge 
of the stellar radius has little effect on the parameter estimation. 

The above analysis implies that with a suitably short exposure time, our formulae 
provide an unbiased estimation of the spin parameters statistically, if we have reasonably 
accurate prior knowledge of the other parameters of the system. In reality, however, it may 
be more relevant to ask about the reliability of the confidence level of the spin parameters 
derived from a single realization, rather than from the 50,000 realizations. This is related to 
the above analysis statistically, but perhaps it is more appropriate to the situation that one 
encounters in any observation. For this purpose, we randomly select one from the 50,000 




(54) 
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realizations and compute 1 and 2 a confidence contours from the relative confidence levels, 

A X 2 = x\VsmI s ,\) -x 2 (\/sinJ Si0 ,Ao), (55) 

where VsinI S fi and Ao are the best-fit values. Figure 13 shows the estimated parameter 
regions on the Vsinl s versus A plane, and the best-fit values are indicated by the plus signs 
in each panel. The corresponding radial velocity curves are depicted in Figure 14, together 
with both the best-fit and the true curves (solid and dashed curves, respectively). 

Basically, Figure 13 demonstrates that the spin parameters V sin I s and A can be con- 
strained around the best-fit values with the 1 a errors of A^sin/ S ~ ±300m s _1 and 
A A ~ ±4°, which greatly improves the uncertainties (by a factor of 4) relative to the previ- 
ous result by Queloz et al. (2000, see Table 3). 



7. Conclusions and discussion 

We have discussed a methodology to estimate the stellar spin angular velocity and its 
direction angle with respect to the planetary orbit for transiting extrasolar planetary systems 
using the RM effect previously known in eclipsing binary stars (Rossiter 1924; McLaughlin 
1924; Kopal 1990). In particular we have derived analytic expressions of the radial velocity 
anomaly, Av s , which are sufficiently accurate for extrasolar planetary systems. If the stellar 
limb darkening is neglected, the expression is exact. We have extended the result to the 
case with limb darkening and obtained approximate but accurate analytic formulae. For a 
typical value of 7 = R p /R s ~ 0.1, the formulae reduce to a simple form (eqs.[40], [43], [44], 
[45], [48], and [49]): 

Av s = n s X p sin I s 7 _* 2 _ l - 2J / (56) 
1 7 e 1 3 7 I 



during the complete transit phase and 



;i- e ){-2 C + 7 2 cos- 1 (C/7)} + r 



Wa 



77 

Av s = Q s X p sinI s —r- n , — - — ~— f x (57) 

during the egress/ingress phases, where 

W2 . ^- X {-^'\ (58) 

Wi = I^ 2 _ ir{l _ ()Xc _0iI!^l_ y (59) 
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where g(x; a, b) is defined in equation (A17). The definition and the meaning of the variables 
in the above expressions are summarized in Table 1. 

The numerical accuracy of the above formulae was checked using a specific example of 
the transiting extrasolar planetary system, HD 209458, and we found that they are accurate 
within a few percent. Our analytic formulae for the radial velocity anomaly are useful in 
several ways. One can estimate the planetary parameters much more efficiently and easily, 
since one does not have to resort to computationally demanding numerical modeling. Fur- 
thermore, the resulting uncertainties of the fitted parameters and their correlations are easily 
evaluated. To illustrate these advantages specifically, we performed a parameter estimation 
applying the formulae against mock data for HD 209458. We showed that with precision 
data obtainable by 8-10 m class telescopes such as Subaru HDS, our formulae improve the 
efficiency and robustness of estimating the spin parameters, Vsinl s and A. Furthermore, 
the combined data analysis with asteroseismology (e.g., Gizon & Solanki 2003) and/or the 
correlation between the mean level of emission and the rotation period (e.g., Noyes et al. 
1985) may break the degeneracy between V and I s in extrasolar planetary systems. 

Among the recently discovered transiting extrasolar planetary systems, i.e., TrES-1 by 
the Trans-Atlantic Exoplanet Survey (Alonso et al. 2004) and OGLE-TR 10, 56, 111, 113, 
132 by the Optically Gravitational Lens Event survey (e.g., Udalski et al. 2002a,b,c, 2003; 
Konacki et al. 2003; Bouchy et al. 2004; Pont et al. 2004), TrES-1 has similar orbital period 
and mass to those of HD 209458b, but its radius is smaller. Thus, it is an interesting target 
to determine the spin parameters via the RM effect; if its planetary orbit and the stellar 
rotation share the same direction as discovered for the HD 209458 system, it would be an 
important confirmation of the current view of planet formation out of the protoplanetary 
disk surrounding the protostar. If not, the result would be more exciting and even challenge 
the standard view, depending on the value of the misalignment angle A. 

We also note that the future satellites Corot and Kepler will detect numerous transiting 
planetary systems, most of which will be important targets for the RM effect in 8 - 10 m 
class ground-based telescopes. We hope that our analytic formulae presented here will be a 
useful template in estimating parameters for those stellar and planetary systems. 

Finally, it is interesting to note that the RM effect may even be used as yet another new 
detection method of transiting planetary systems. For the HD 209458 system, the stellar 
radial velocity amplitude due to the Kepler motion is around 100 m s -1 . Since the stellar 
rotation velocity is around 4 km s _1 , the maximum radial velocity anomaly due to the RM 
effect is ~ 40(7/0. l) 2 m s _1 and thus is very close to the former. On the other hand, the 
latter could be significantly larger if the host star rotates faster, and/or the host star (the 
planet) has a smaller (larger) radius. In extreme cases, the radial velocity curve, with a 
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periodicity of a few days for instance, is barely detectable, but the velocity anomaly is quite 
visible for a few hours of the transiting phase. Thus, the conventional radial velocity curve 
analysis might have missed some of the potentially interesting spectroscopic signature of 
transiting planets. A search for such specific signatures deserves an attempt against existing 
or future spectroscopic samples that do not show any clear conventional feature of radial 
velocity periodicity. 

In conclusion, we have demonstrated that the radial velocity anomaly due to the RM 
effect provides a reliable estimation of spin parameters. Combining data with the analytic 
formulae for radial velocity shift Av s , this methodology becomes a powerful tool in extracting 
information on the formation and the evolution of extrasolar planetary systems, especially 
the origin of their angular momentum. Although it is unlikely, we may even speculate that 
a future RM observation may discover an extrasolar planetary system in which the stellar 
spin and the planetary orbital axes are antiparallel or orthogonal. This would have a great 
impact on the planetary formation scenario, which would have to invoke an additional effect 
from possible other planets in the system during the migration or the capture of a free- 
floating planet. While it is premature to discuss such extreme possibilities at this point, 
the observational exploration of transiting systems using the RM effect is one of the most 
important probes for a better understanding of the origin of extrasolar planets. 

We thank Christopher Leigh, Norio Narita, and Edwin Turner for useful discussions 
and comments. We are also grateful to Didier Queloz for providing the ELODIE data of the 
radial velocity of HD 209458. This work is supported in part by a Grant-in- Aid for Scientific 
Research from the Japan Society for Promotion of Science (14102004, 16340053). 
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A. Approximation of integrals 

In this appendix, the approximate expressions for the integrals Wi (i — 1 ~ 4) defined 
in section 5 are derived. Below, we present the results separately for the integrals W\ and 
W2 in Appendix A.l and for W3 and W4 in Appendix A. 2. 



A.l. Integrals W 1 and W 2 



First, consider the integral W±. For further reduction of the integral, it is convenient to 
rewrite equation (41) in terms of the new variables 

(x,z) = i? p <r(cos <p, sin p) + (X P ,Z P ), (X P ,Z P ) = R s p(cos9 : sm9). (Al) 

Then we obtain 



W 1 = - [ da [ dip a\[l - p 2 sin 2 (0 - <p) - {70- + pcos (0 - p)} s 
n Jo Jo v 



(A2) 



Here, the variable p runs from to 1 — 7. In the above expression, the integral over a is 
analytically expressed, and a part of it is further integrated out. The resulting expression 
becomes 



W 1 



2 



:i-p 2 ) 1/2 -va + wb, 



(A3) 



where wa and wb are given by 

o2tT 



w A = 



7T 7 



1 f 

- d(p(l - p 2 - 7 2 - 2p7 cos<p) 3/2 
3 Jo 

1 /" 27r 

+ - / d<p (p 2 cos (/? + P7) cos <p(l — p 2 — 7 2 — 2p7 cos <p) 

2 Jo 



1/2 



,(A4) 



1 /" 27r 

= - — x / d<p pcos <p(l - p 2 sin 2 p) 

2tt7 2 y 



x < sin 



-1 



pcosp 



a/1 — p 2 sin 2 p 



_i / 7 + pcosp 
sin ' 



a/1 — p 2 sin 2 p 



■(AS) 



Note that the integrals u>a and »b are analytically expressed in terms of the complete elliptic 
integrals of the first kind K(m), the second kind E(m), and the third kind H"(n, m): 

- = ^<i-h + P) 2 }" 2 

x [{3(1 - p>Y - 7 2 (71 - 17/) + 68 7 4 } B(m) 



- 21 - 



w B 



2 



-{(7 - P? ~ l}(3p 2 - 3 + 8 7 2 ) K{m)\ , (A6) 
1-(1-P 2 ) 1/2 (l-0 



2 

1 

4 



45tt7 4 [1 - (7 + p) 2 ] 1 ^ 

x { [2 7 6 + (29 - 7p 2 ) 7 4 + (26 - 34p 2 + 8p 4 ) 7 2 + 3(1 - p 2 ) 3 ] X(m) 

+ [1 - (p + 7 ) 2 ] [2 7 4 + (31 - 7p 2 ) 7 2 - 3(1 - p 2 ) 2 ] E(m) } 
2 



3vr 7 2 [l - ( 7 + p) 2 ]^ 



{[1 - n(p + 7 )/2p] U(n, m) + c.c.} , (A7) 



where c.c. denotes the complex conjugate, n = 2p [p — i(l — p 2 ) 1 / 2 ], and m = — 4 7 p/[l — 
( 7 + p) 2 ]. To evaluate equations (A6) and (A7), careful treatments are required at the edge 
p = 1 — 7 , where the argument of the elliptic integral and coefficients of K{m) and H(n, m) 
in % apparently diverge. Since we are concerned with planetary systems with a small ratio 
R p <^ R s , it is practically useful to derive the approximate expressions. In this case, we do 
not have to use the complicated expressions (A6) and (A7), but can expand the integrands 
(A2) and/or (A4), (A5) in powers of 7 = R p /R s . Then each term in the expansion with 
respect to 7 can be analytically integrated. The results up to the second order in 7 become 

2 M 21/2 ^ 1 2 ^ 1 15p 4 -28p 2 + 16 



WA * W {l ~ P) VS) - 16 (1-P 2 )3/ 2 

- 7 2 J- 17p 6 -64p 4 + 104p 2 -32 

7 128 (l-p 2 ) 7 /2 + u ^ >> ^ 

4 (1-p 2 ) 3 / 2 7 128 (l-p 2 ) 7 /2 + u ^ >■ ^> 

Substituting these expressions into equation (A3) and collecting the terms in powers of 7, 
we finally obtain 

Wi - (1 - P 2 ) 1/2 - 7 2 8{ *:^ )3/2 + 0{l% (A10) 

Next, consider the integral W 2) whose analytic expression is also obtained through the 
same procedure as mentioned above. Using equations (Al), one writes equation (42) as 

W 2 = W 1 + w c , (All) 

where wc is given by 

wc — - — 77 da dip a 2 coscpy 1 — p 2 sin 2 (9 — <p) — {70- + pcos (9 — y?)} 2 . (A12) 

Tip COS J J Q V 
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In the above expression, the integral over a is analytically performed, and the resulting 
expression for the integrand is expanded in powers of 7. Further integrating it over ip, we 
obtain 

^-7 2 4(1 _ p2)1/2 +<?(7 4 ). (A13) 
Thus, the perturbative expansion for Wi becomes 

P 2 ) 1/2 - 7 2 8( i Z ^3/2 + °^)- ( Al4 ) 

In Figure 15, to check the validity of the perturbation results, the approximate ex- 
pressions for the integrals W\ and W 2 are plotted as a function of p. The results are then 
compared with those obtained from the numerical integration. As is expected, the pertur- 
bative expressions (A10) and (A14) give a quite accurate approximation as long as the ratio 
of the planetary radius to the stellar radius is small. Note that the approximation is even 
better for the slightly larger value 7 = 0.3, with a few percent level for the fractional error. 



A. 2. Integrals W 3 and W A 

As for the integrals W3 and W 4 given by equations (46) and (47), one can partially 
evaluate the integrals with the knowledge of the integral region (32). The resulting forms 
are summarized as 



7T f X ° 

W 3 = -(l-x ) 2 (2 + x )+ dxg{x\ri P ,i), (A15) 

7r r xo 

W 4 = -(l-x 2 ) 2 + dxx g(x; Vp , 7), (A16) 

5 Jxo+C-l 

where we defined the function g(x; i] p , 7) as 

t~ \ n ■ -1 [ 7 2 -(^-l-?7 P ) 2 l 1/2 
9{x;V P ,V = (l-ar)sin <^ ^— ^ > 

+ ^2(1 + np ){x - x){^ - (x - 1 - Vp ) 2 }. (A17) 

Since the above one-dimensional integrals cannot be evaluated analytically, one may derive 
an approximate expression applicable to the 7 < 1 cases. Note, however, that a naive 
treatment by the perturbative expansion regarding 7 as a small parameter can break down 
at rj p ~ —7. Even in the r\ p > 7 case, perturbative expression gives a worse approximation. 
For an accurate evaluation of the integrals, a more dedicated treatment other than the 
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perturbative expansion is required. One clever approach, which we adopt here, is to replace 
the function g with 



, (a — x)(x — b) /, 10 \ 

9{X;V P ,V ► 9{X C ; V P , V \ 1 w JT, (A18) 

" (a - x c )(x c - b) 

where we set a = xq + ( — 7, b = xq, and x c = (a + b)/2, and to use the integral formula 



J a 



i (a — x)(x — b) 7T /T , /,-,„x 

^/ >LJL->) = 4 (t - a) - (Al9) 



Then, one can approximate the one-dimensional integrals in equations (A15) and (A16) as 

7T 

dx g(x; r? p , 7) — > - (7 - C) #(x c ; ?? p , 7), (A20) 
-7 4 

7T 

dxxg(x;r] p , 7) — ► - (7 - C)(2x + C - 7) 9( x c, %, l)- (A21) 



As long as the planetary transit system has 7 <C 1, the above expression in fact gives an 
accurate prescription. Note, however, that a naive use of the formulae (A20) and (A21) 
leads to inconsistent radial velocity curves, which do not satisfy the junction condition at 

V P = -r- 

W 3 = tt 1 2 W^| p=1 _ 7 , W 4 = vr 7 2 (l - 7) W 2 \ p=1 ^. (A22) 

To preserve the consistency, we modify equations (A20) and (A21) so as to satisfy the 
junction condition (A22) by multiplying the numerical factor. The final expressions for W3 
and W4 become 



W 3 ~ |(1 - x,f{2 + s ) + I 7(7 - C) jt ^ - 7) , ~ 7), (A23) 

6 2 5(1-7; -7, 7) 



W 4 ~ ^(1 - x ) 2 (l + x ) 2 + f 7(7 ~ ^ ,f Xc] Vp ' 7) , W 2 (l ~ 7). (A24) 
8 2 g(l -7; -7, 7) 

In figure 16, substituting equations (A10) and (A14) into the above results, the approx- 
imate expressions for W 3 and W 4 are depicted as a function of i] p . Although it is a tricky 
treatment, it turns out that the expressions (A23) and (A24) give an accurate approximation 
for the range of our interest. Compared to the integrals W\ and W 2 , the fractional error of 
the approximations seem slightly large; however, the contribution of the integrals W3 and 
W4 to the radial velocity shift Av s is relatively small. With a typical parameter 7 ~ 0.1, the 
resulting fractional error remains only a few percent. 
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Table 1. List of notation 



Variables Definition Meaning 



Orbital Parameters 



m p 


Sec.2 


Planet mass 


m s 


Sec.2 


Stellar mass 


a 


Fig.l 


Semimajor axis 


e 


Fig.l 


Eccentricity of planetary orbit 


w 


Fig.l 


Negative longitude of the line of sight 


i 


Fig.2 


Inclination between normal direction of orbital plane and y-axis 




Eq.[l] 


Distance between star and planet (see Fig.l) 


f 


Eq.[2] 


True anomaly (see Fig.l) 


E 


Eq.[2] 


Eccentric anomaly 


n 


Eq.[3] 


Mean motion 


M 


Eq.[4] 


Mean anomaly 


Internal Parameters of Star and Planet 


/. 


Fig.2 


Inclination between stellar spin axis and y-axis 


A 


Fig.3 


Angle between z-axis and normal vector h p on (x, z)-plane 


Ms 


Eq.[12] 


Angular velocity of star (see Fig.2) 


R s 


Sec.4 


Stellar radius 


R p 


Sec.4 


Planet radius 


e 


Eq.[38] 


Limb darkening parameter 


V 


Sec. 6 


Stellar surface velocity, R S Q S 


Mathematical Notation 


X p 


Sec.4 


Position of the planet 


7 


Eq.[25] 


Ratio of planet radius to stellar radius, R p / R s 


V P 


Eq.[28] 


See Fig.6 


x 


Eq.[33] 


See Fig.6 
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Table 1 — Continued 



Variables 


Definition 


Meaning 




Eq.[33] 


See Fig.6 


C 


Eq.[34] 


See Fig.6 



Table 2. Parameters of the HD 209458 system 



Parameters 


Estimated Values 


Fractional Errors (%) 


m s 


1.1 ±0.1 M a 


9.1 


m p 


0.69 ±0.05 M/ 


7.2 


e 


0, b - 0.1 c 




a 


0.0468 a ±0.0014 d AU 


3.0 


zu 


100° c 




Porb 


3.52474 ± 0.00007 days a 


0.002 


R s 


1.146 ±0.050 R a 


4.4 


R p 


1.347 ±0.060 R/ 


4.5 


i 


86°. 1 ±0°.l b 


0.1 


e 


0.64±0.03 a 


4.7 



a From Brown et al. (2001). 

b From http:/ /www. obspm.fr/encycl/HD209458. html. 
c From http://exoplanets.org/esp/hd209458/hd209458.html. 
d This error is calculated from those of m s , m p , and P or b- 
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Table 3. Spin Parameters Derived by Queloz et al. (2000) and Notational Differences 



This Paper Queloz et al. (2000) Best-Fit Solutions a 



V sin I s v sin i 


3.75 ± f. 25 b 


sgn(A) cos 1 {cosAsini} 


a 


±3.9+2i C 


A sgn(a) cos -1 


cos a 


for (a, tt p ) = (3.9, 0) c 


(l— cos 2 Up sin 2 a) ^ 




21.7 for (a,Q p ) = (22, 100) c 
-24.7 for (a, = (-25,80) c 



a The parameter ft p is related to the inclination angle and is constrained through 
cos ftp = — cos(86?l)/ cos a. 

b In units of km s _1 . 

c Units for a, A, ftp are degrees. 
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observer's 
direction 



Fig. 1. — Schematic top view of the planetary orbit. The star is located at the focus point 
F (see Table 1 for meaning of symbols). 
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z 




Fig. 2. — Schematic configuration of the stellar spin axis, the planetary orbital plane, and 
the observer's line of sight (see Table 1 for meaning of symbols). 




Fig. 3. — Projected view of the planetary orbital plane from the line of sight (F-axis in this 
case) . 
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time 



Fig. 4. — Schematic illustration of the velocity curve anomaly due to the RM effect. The 
four different paths of a planet, a-d, in the left panel correspond to the velocity curves in 
the right panel. 




Fig. 5. — Planetary transit: ingress, complete transit, and egress phases. 



-32 - 




Fig. 6. — Schematic illustration of the configuration of the system at ingress of the planet 
in the new coordinates. 
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Fig. 7. — Analytic radial velocity curves for the RM effect without stellar limb darkening 
(e = 0). 
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Fig. 8. — Same as Fig. 7, but with a linear limb-darkening effect (e = 0.64). 
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Fig. 9. — Variation of radial velocity shift with respect to a parameter variation p — > p+dp as 
a function of time. The vertical axis denotes the variation of radial velocity shift normalized 
by the fractional error in each parameter dp/p ( eqs. [52] and [53] for left and right panels, 
respectively). These quantities are evaluated around the fiducial values for the HD 209458 
system summarized in Table 2 using the analytic formulae for the RM effect. 
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Fig. 10. — Comparison of different models of the radial velocity anomaly curves for the HD 
209458 system without limb darkening (e = 0). The top panel shows the radial velocity shifts 
obtained from the numerical evaluation of expression (20) (solid curve) and the analytic for- 
mulae with and without the effect of exposure time (dotted and dashed curves, respectively) . 
The bottom panel shows the residuals between the two radial velocity shifts among three 
curves. 
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Fig. 11. — Same as Fig.10, but with limb darkening (e = 0.64). 
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Fig. 12. — Joint probability distribution of the best-fit parameters (Vsin J s , A) in the like- 
lihood analysis for the mock HD 209458 observation. We adopt the prior assumption for 
fiducial parameters (top left), R p = 1.407 Rj (top right), R s = 1.196 R Q (bottom left) and 
a = 0.0482 AU (bottom right). In each panel, the cross indicates the correct values of 
(Vsml s , A), and the contour levels around it represent the 1 and 2 a levels of the probabil- 
ity distribution. The solid lines projected onto the horizontal and vertical axes represent the 
probability distributions of V sin I s and A, respectively. Note that all the mock data created 
in the likelihood analysis assume the parameters listed in Table 2. Thus, except for the top 
left panel, the difference between the incorrect prior assumption and the correct value for 
(R p , R s , a) was used in each panel. 
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Fig. 13. — Expected constraints on the spin parameters (Vsinl s , A) from the x 2 fit to 
the mock data of radial velocity curves for the HD 209458 system. The thin and the thick 
lines show the results obtained from the different mock samples created according to the 
same parameters, as listed in table 2. The two contour curves in each line width represent 
the 1 and 2 a confidence levels. The locations of the crosses represent the true values of 
the spin parameters (Vsin/ S , A), while the locations of the plus signs indicate the best- 
fit parameters, which were estimated under the prior assumption for the limb-darkening 
coefficient and the orbital eccentricity as the corect parameters (top left), AR p = 0.06 Rj 
(top right), AR S = 0.05 R Q (bottom left), and Aa = 1.4 x 10~ 3 AU (bottom right). 
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Fig. 14. — Comparison between the mock data for the radial velocity curves used in Fig. 13 
and the best-fit curves for the parameters (Vsml s , A) (solid curves) determined under the 
prior assumption of (e, e) for the fiducial parameters (top left), AR P = 0.06 Rj (top right), 
AR S = 0.05 = Rq (bottom left), and Aa = 1.4 x 10~ 3 AU (lower-right). The result depicted 
in the upper and lower window of each panel corresponds to the thin and thick contour lines 
in Fig. 13, respectively. For comparison, we also plot the correct radial velocity curves in 
each plot (dashed curves). 
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Fig. 15. — Top windows: Approximate vs. numerical evaluations of the integrals W\ {left) 
and W2 (right) as a function of p = [X p + ZTj / R s . The solid lines represent the numerical 
evaluation of equations (A3) and (All), while the dashed lines are the approximation based 
on the perturbative expansions (A10) and (A14). Note that the variable p runs from to 
1 — 7. Bottom windows: Fractional errors W^ approx ^ /W^ num ^ — 1 for the integrals W\ (left) 
and W 2 (right). 
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Fig. 16. — Top windows: Approximate vs. numerical evaluations of the integrals W3 (left) 
and W4 (right) as a function of r) p . The solid lines represent the numerical evaluation of 
equations (A15) and (A16), while the dashed lines are the approximation based on the 
perturbative expansions (A23) and (A24). Note that the variable rj p runs from —7 to +7. 
Bottom windows: fractional errors 

^(approx) / W (num) _ 1 for the integra l s ^ W 4 

(right). 



